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Dissociation products and structures of solid H2S at strong compression 
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Hydrogen sulfides have recently received a great deal of interest dne to the record high supercon¬ 
ducting temperatures of up to 203 K observed on strong compression of dihydrogen sulfide (H 2 S). 

A joint theoretical and experimental study is presented in which decomposition products and struc¬ 
tures of compressed H 2 S are characterized, and their superconducting properties are calculated. In 
addition to the experimentally known H2S and H3S phases, our first-principles structure searches 
have identified several energetically competitive stoichiometries that have not been reported pre¬ 
viously; H2S3, H3S2, and H4S3. In particular, H4S3 is predicted to be thermodynamically stable 
within a large pressure range of 25-113 GPa. High-pressure room-temperature X-ray diffraction 
measurements confirm the presence of H3S and H4S3 through decomposition of H2S that emerge at 
27 GPa and coexist with residual H2S, at least up to the highest pressure studied in our experiments 
of 140 GPa. Electron-phonon coupling calculations show that H4S3 has a small Tc of below 2 K, 
and that H 2 S is mainly responsible for the observed superconductivity of samples prepared at low 
temperature (<100K). 

PACS numbers: 62 . 50 .-p, 61 . 50 .Ah, 61 . 05 .cp, 74 . 70 .Ad 


Superconductivity with a transition temperature 
of up to 203 K was observed recently in solid H2S at 
Megabar pressures, which is the highest record among 
all known superconductors [1]. Ashcroft suggested that 
metallic hydrogen would be a superconductor at high 
pressures with a Tc around room temperature [2], and 
subsequently predicted that hydrogen-rich metallic com¬ 
pounds might also be superconducting at high pressures 
[3]. Early theoretical studies focussed on high-pressure 
silicon and aluminum hydrides [4, 5] , and a number of po¬ 
tential high-temperature superconductors have now been 
proposed in compressed hydrogen-rich compounds, with 
TcS estimated in the range 40-250 K (e.g., GaHa [6, 7], 
SnH4 [8], GeH4 [9], NbH4 [10], SiaHg [11], SiH4(H2)2 [12], 
GaHe [13], YH3 [14], YH4 and YHg [15]). They have not 
been realized in practice, in part because of demanding 
experimental challenges. 

The high pressure phase diagram of H2S has been stud¬ 
ied extensively. H2S is a sister molecule of H2O, and 
is the only known stable compound in the H-S system 
at ambient pressure. High-pressure diamond anvil cell 
experiments led to the discovery of a metallic phase at 
about 96 GPa [16-27]. However, partial decomposition of 
H2S and elemental sulfur was observed in Raman [26] and 


XRD studies [21] at room temperature above 27 GPa. 
H2S had not been considered as a candidate for super¬ 
conductivity because it was believed to dissociate into 
elemental sulfur and hydrogen under high pressures [28] . 
Recent first-principles structure searches predicted ener¬ 
getically stable metallic structures of H2S above 110 GPa 

[29] and excluded dissociation into its elements. An esti¬ 
mated maximum Tc for metallic H2S of 80 K at 160 GPa 
was predicted [29]. Motivated by this study, Drozdov et 
al. [1] performed high-pressure experiments on solid H2S 
looking for superconductivity and found an astonishing 
Tc of 203 K at 155 GPa [1]. H2S shows complex supercon¬ 
ducting behavior at high pressures with the emergence of 
two different superconducting states. Samples prepared 
at low temperature (100 K) have a Tc of ^30 K at 110 
GPa at the onset of superconductivity, which increases 
rapidly to a maximum value of 150 K at 200 GPa, while 
samples at room temperature or above show a maximum 
Tc of 203 K at 155 GPa. 

Strobel et al. synthesized another H-S compound, H3S, 
by compressing a mixture of H2S and H2 above 3.5 GPa 

[30] . The superconducting Tc of H3S at 200 GPa was 
recently predicted to be as high as 191-204 K [31], with 
H3S in a cubic Im3m structure, which is known already 
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in H 3 O at terapascal pressures [32] through decomposi¬ 
tion of H 2 O. The agreement between experimental [1] and 
theoretical values of [31] led to the proposal by Droz¬ 
dov et al. that H 3 S could be formed by decomposition of 
H 2 S and might be responsible for the observed supercon¬ 
ducting state at 203 K [1] . This proposal was supported 
by first-principles density-functional-theory (DFT) stud¬ 
ies which suggested that it is thermodynamically favor¬ 
able for H 2 S to decompose into H 3 S -I- S [33-35] at pres¬ 
sures above 43 GPa [34]. Apparently, there is an urgent 
need to characterize the decomposition products of com¬ 
pressed H 2 S in an effort to build an understanding of the 
complex superconducting behavior exhibited by the H-S 
system. 

Here we present a joint theoretical and experimen¬ 
tal study of compressed H 2 S which clarifies the possi¬ 
ble decomposition products and their structures. First- 
principles DFT structure searches were used to predict 
several new stoichiometries (H 2 S 3 , H 3 S 2 , and H 4 S 3 ) and 
a new structure in H 3 S not reported previously. Room- 
temperature high pressure X-ray diffraction (XRD) ex¬ 
periments demonstrate that above 27 GPa, H 2 S partially 
decomposes into S -I- H 3 S -I- H 4 S 3 . H 4 S 3 emerges as the 
major component at around 66 GPa and coexists with a 
small fraction of H 3 S and residual H 2 S, at least up to the 
highest pressure studied experimentally of 140 GPa. 

Extensive structure searches over 44 H-S stoichiome¬ 
tries at 25, 50, 100 and 150 GPa were performed us¬ 
ing the GALYPSO [36, 37] and AIRSS [4, 38] methods, 
which have been successfully used to investigate struc¬ 
tures of materials at high pressures [4, 13, 39-45]. The 
underlying structural relaxations were performed using 
the Vienna ab initio simulation package (VASP) [46] 
for CALYPSO and the CASTEP plane-wave code [47] 
for AIRSS. Electron-phonon coupling (EPC) calculations 
were performed with density functional perturbation the¬ 
ory using the Quantum-ESPRESSO package[48]. XRD 
data were collected at the 15U1 beamline at the Shanghai 
Synchrotron Radiation Eacility (SSRF) with a monochro¬ 
matic beam of wavelength 0.6199 A. The diffraction pat¬ 
terns were integrated with the E1T2D computer code 
[49] and htted by Rietveld profile matching using the 
GSAS-l-EXPGUl programs [50, 51]. More information 
about the calculations and experiments is provided in 
the Supplemental Material [52]. 

Figure 1 shows convex hull diagrams at 25, 50, 100 and 
150 GPa which summarize the results of the structure 
searches. The effects of including quantum zero-point 
vibrational motion are significant and they tend to in¬ 
crease with pressure. Our results suggest that up to five 
H-S compounds lie on the convex hull at some pressures 
and are therefore thermodynamically stable. Besides the 
experimentally known H2S and H3S compounds, we pre¬ 
dict three additional stable compounds: H4S3, H3S2 and 
H2S3. Note that H2S is theoretically found to be sta¬ 
ble only below 25 GPa, while H3S appears at all pres- 
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FIG. 1: (Color online) Results from structure searching at 
25 (a), 50 (b), 100 (c) and 150 GPa (d). Gonvex hulls are 
shown as continuous lines, with (red) and without (black) the 
inclusion of zero point vibrational enthalpy (ZPE). 


sures considered. Enthalpy calculations show that H 3 S 
becomes energetically more stable than H 2 S-I-I/ 2 H 2 at 
around 6 GPa. The newly predicted H 3 S 2 and H 2 S 3 
phases have very narrow pressure ranges of stability and 
are unstable above 34 and 25 GPa, respectively (Supple¬ 
mental Material, Fig. SI [52]). We have therefore omitted 
further discussion of these compounds. The correspond¬ 
ing crystallographic parameters and phonon dispersion 
curves are provided in the Supplemental Material [52]. 

For H 3 S, besides the PI, Gccm, R3m and lm3m struc¬ 
tures of earlier studies [31], our searches predict a mon¬ 
oclinic C2/c structure (4 f.u/cell) that was not reported 
earlier. The C2/c structure consists of weakly bonded 
H 2 S and H 2 molecules (Supplemental Material, Fig. S7b 
[52]), and is calculated to be more stable at pressures of 
2-112 GPa than the Gccm structure proposed previously 
[31] (Supplemental Material, Fig. S7a [52]). Static-lattice 
enthalpy calculations give a zero-temperature phase se¬ 
quence for H 3 S of PI ^ C2/c (2 GPa) R3m (112 GPa) 
^ Im3m (180 GPa). 

H 4 S 3 adopts an orthorhombic P2i2i2i (4 f.u./cell) 
structure that consists of weakly bonded HS and H 2 S 
molecules at 25 GPa (Fig. 2a). The H-S bond lengths 
within the HS and H 2 S molecules are 1.354 A and 1.387- 
1.391 A, respectively, which are significantly shorter 
than the H-S separation of 1.913-1.932 A between 
molecules. With increasing pressure, the neighboring 
molecules bond with each other forming planar H-S-H- 
S zigzag chains and puckered H-S-H-S chains, respec¬ 
tively. P2i2i2i transforms to a Pnma structure at 60 
GPa (Fig. 2c). The convex hull data suggests two syn¬ 
thesis routes for H 4 S 3 : (i) decomposition of 8 H 2 S—^-S + 
4 H 3 S -I- H 4 S 3 above 25 GPa; (ii) reaction of 4 H 3 S -I- 5S 
—>■ 3 H 4 S 3 in the pressure range of 25-113 GPa. Theo- 
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FIG. 2: (Color online) Energetically favorable structures of 
H4S3 with space groups P2i2i2i (a) and Pnma (b). Small 
and large spheres represent H and S atoms, respectively, (c) 
Calculated enthalpy curves of the P2i2i2i structure with re¬ 
spect to the Pnma structure for H4S3 as a function of pressure, 
(d) Decomposition enthalpy curves of H4S3 into (H3S-I-S) as 
a function of pressure. 


retically, it is found that H 4 S 3 decomposes into H 3 S -b S 
above 113 GPa (Fig. Id and Fig. 2d). 

Powder XRD patterns obtained on increasing pressure 
from 10 to 140 GPa at room temperature are shown 
in Fig. 3a. The XRD patterns collected at pressures 
up to 46 GPa are in excellent agreement with previous 
data [17-21], and the successive transitions of phase F —)■ 
phase IV ^ phase V are well reproduced. The XRD data 
at 10 GPa correspond to phase I'. Phase IV with addi¬ 
tional peaks at around 12° and 15° (shown by arrows in 
Fig. 3a) was observed at 16 GPa. The diffraction peaks 
of phase IV weaken at pressures above 27 GPa, and a 
new diffraction profile observed at 46 GPa is responsible 
for phase V. The XRD data show that the IV^-V tran¬ 
sition begins above 27 GPa, in excellent agreement with 
previous results [18, 21 ]. 

Previous high-pressure Raman [26] and XRD [21] stud¬ 
ies have claimed that decomposition of H 2 S occurs at 
room temperature above 27 GPa. Indeed, we found that 
H 2 S partially decomposes in phase V. Unfortunately, we 
have not found it possible to resolve the decomposition 
products and their crystal structures from our current 
XRD data. Therefore we use the predicted structures 
and convex hull data to help in analysing the experimen¬ 
tal data. At 50 GPa, our calculations suggest an ener¬ 
getically allowed dissociation path of 8 H 2 S —+ 4 H 3 S + 
H 4 S 3 (Fig. lb). The XRD profile at 46 GPa was there¬ 
fore fitted to a mixture of H 2 S -b S -b H 3 S -b H 4 S 3 by 
performing Rietveld refinements using the most energet¬ 
ically stable structures. Remarkably, we found that the 
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FIG. 3: (Golor online) (a) XRD patterns of H 2 S collected at 
various pressures and at room temperature with an incident 
wavelength of 0.6199 A. [(b) and (c)] Rietveld refinements of 
XRD profiles at 45.6 GPa based on the compositions PC-H 2 S, 
I4i/acd-S, G 2 /C-H 3 S and P 2 i 2 i 2 i-H 4 S 3 with phase fractions 
of 1147:85:31:1, and at 65.7 GPa based on compositions of 
Pnma-Il 4 S 3 , G 2 /C-H 3 S and PC-H 2 S with phase fractions of 
5:3.4:1, respectively. The cross symbols and red solid lines 
represent observed and fitted patterns, respectively. The solid 
lines at the bottom are the difference between the observed 
and fitted patterns. Vertical bars under the pattern repre¬ 
sent the calculated positions of reflections arising from the 
compositions. 


XRD profile can be well indexed by a mixture of Pc- 
structured H 2 S, I4i/acd-structured S, G2/c-structured 
H 3 S and P2i2i2i-structured H 4 S 3 , with phase fractions 
(ratios of numbers of unit cells) of 1147:85:31:1, as shown 
in Fig. 3b. The existence of a large proportion of H 2 S 
demonstrates the partial decomposition. We also at¬ 
tempted other Rietveld refinements fitting the XRD pat¬ 
terns to pure H 2 S, H 2 S -b S -b H 3 S or S -b 4 H 3 S -b H 4 S 3 , 
but all of these fits gave poorer results (Supplemental 
Material, Fig. SIO [46]). The calculated decomposition 
pressure (30 GPa) for 8 H 2 S —> S -b 4 H 3 S -b II 4 S 3 (Supple¬ 
mental Material, Fig. S13 [52]) is in excellent agreement 
with the value of 27 GPa observed in experiment [21]. 

With increasing pressure, more H 2 S decomposes and 
its contribution to the XRD signal is reduced. The XRD 
pattern collected at 66 GPa shows entirely different fea¬ 
tures to that at 46 GPa. Rietveld refinement shows that a 
mixture of Pnma-structured H 4 S 3 , C2/c-structured H 3 S, 
and Pc-structured H 2 S with phase fractions 5:3.4:1 gives 
the best fit to the experimental data (Fig. 3c). The dis¬ 
appearance of elemental S and the reduction in the ratio 
of H 3 S are understandable since a reaction of 4 H 3 S -b 
5S —> 3 II 4 S 3 takes place as inferred from our convex hull 
calculations (Fig. lb and Ic). The two shoulders on the 
main peak at 15° originate primarily from H 3 S and H 2 S. 
These shoulders weaken when the pressure is increased 
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FIG. 4: (Color online) Rietveld refinements of XRD profiles 
at 81.3 GPa based on Pnma-H4S3 + C 2 /C-H 3 S + Pmc2i-H2S 
with phase fractions of 43:6:1 (a) and at 139.5 GPa based on 
Pnma-H 4 S 3 + R3m-H3S + P-I-H 2 S with phase fractions of 
56:7:1 (c). 


to 82 GPa and the XRD pattern can then be well in¬ 
dexed by Pnma-structured H 4 S 3 , C2/c-structured H 3 S 
and Pmc2i-structured H 2 S with a major contribution 
( 86 %) from H 4 S 3 (Fig. 4a). Our results demonstrate that 
H 4 S 3 coexists with H 3 S and H 2 S at least up to 140 GPa, 
the highest pressure studied experimentally (Fig. 4b). 
At this pressure, a refinement based on Pnma-structured 
H 4 S 3 and G2/c-structured H 3 S leads to poorer fits with 
higher Rp and Ru,p values (Supplemental Material, Fig. 
S12 [52]), which supports the existence of residual H 2 S 
(about 1 . 6 %). 

H 4 S 3 becomes metallic at 102 GPa (Supplemental Ma¬ 
terial, Fig. S14 [52]). However, the calculated electron- 
phonon-coupling parameter (A = 0.42) is very small at 
140 GPa due to the low density of states at the Fermi 
level of 0.09 eV“^/atom. As a result, the Tc estimated 
from the Allen and Dynes modified McMillan equation 

[53] with /i* of 0.16-0.13 is only 0.75-2.1 K at 140 GPa. 

In Fig. 5, we compare the calculated Tc values for 

H 4 S 3 , H 2 S and H 3 S to experimental data measured in 
compressed H 2 S [ 1 ], where Tc obtained for samples pre¬ 
pared at low and high temperatures are denoted by L- Tc 
and H-Tc, respectively. On the one hand, the observed 
L-Tc [1] at pressures below 160 GPa can only be quan¬ 
titatively reproduced by H 2 S, while the nature of the 
rapidly increasing L-Tc above 160 GPa remains unclear 
because the calculated Tc values of H 3 S, H 4 S 3 , and H 2 S 
are too high, too low, and tending to decrease, respec¬ 
tively. On the other hand, although the values of Tc 
for Im3m structured H 3 S calculated within the harmonic 
approximation are much higher than the observed H- Tc 

[54] , the inclusion of anharmonic effects reproduces the 
H-Tc data above 180 GPa [35] quite well. However, at 
low pressures around 150 GPa, Tc values [54] estimated 
for R 3 m-H 3 S within the harmonic approximation are ^45 
K lower than the observed H-Tc. Meanwhile, the pre¬ 
dicted Tc for R 3 m-H 3 S increases with pressure, in stark 
contrast to the experimental observation of a decrease in 



FIG. 5: (Golor online) Superconducting transition temper¬ 
atures {Tc) calculated for various H-S compounds and ex¬ 
perimental values for compressed H 2 S [1]. Solid squares and 
circles show experimental data from Fig. 1(b) and Fig. 2 (b) 
of Ref. 1, respectively, where different runs are represented by 
the same symbols. The open stars, squares and circles show 
calculated data from Ref. 29, Ref. 35 and Ref. 54, respec¬ 
tively. Open triangles denote calculated Tc for Pnma-H 4 S 3 
in the present study. Note that the Tc values (open squares) 
for ImSm-HsS taken from Ref. 35 include anharmonic effects, 
while all other estimated TcS are calculated within the har¬ 
monic approximation. 


H- Tc- Apparently, further study is greatly needed to dis¬ 
close the steep Tc increase of L- Tc above 160 GPa and 
the high H-Tc at around 150 GPa. 

We find that kinetics plays an important role in de¬ 
termining the experimentally observed H-S structures. 
Theory suggests that H 2 S and H 4 S 3 decompose above 25 
and 113 GPa (Supplemental Material Fig. S13 and Fig. 
2d), respectively. However, H 2 S and H 4 S 3 are observed 
to persist up to at least 140 GPa. Large kinetic barriers 
must therefore play a major role in suppressing decom¬ 
position at high pressures, as has been found in other 
materials [55-58]. 

In summary, through first-principles structure search¬ 
ing calculations, we predict three new stable H-S com¬ 
pounds with stoichiometries H 3 S 2 , H 2 S 3 , and H 4 S 3 and 
a new C2/c structure of H 3 S, enriching the phase dia¬ 
gram of H-S systems at high pressures. The formation 
of H 4 S 3 and H 3 S was confirmed by XRD experiments to 
occur through decomposition of compressed H 2 S above 
27 GPa resulting in the products S -I- H 3 S -I- H 4 S 3 . H 4 S 3 
becomes a major component at around 66 GPa and is 
stable up to at least 140 GPa, with a small fraction of 
H 3 S and residual H 2 S. We have also examined poten¬ 
tial superconductivity of metallic H 4 S 3 via explicit cal¬ 
culations of electron-phonon coupling parameters and the 
superconducting Tc- Our work suggests that kinetically 
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protected H 2 S in samples prepared at low temperature 
is responsible for the observed superconductivity below 
160 GPa. 
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